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Abstract The long wavelength emission of protostellar objects is commonly attributed to a disk of gas and dust 
around the central protostar. In the first stages of disk accretion or in the case of high mass protostars, the disk 
mass is likely to be sufficiently large, so that the disk self-gravity may have an impact on the dynamics and the 
emission properties of the disk. In this paper we describe the spectral energy distribution (SED) produced by a 
simple, non-flaring, self-gravitating accretion disk model. Self-gravity is included in the calculation of the rotation 
curve of the disk and in the energy balance equation, as a term of effective heating related to Jeans instability. 
In order to quantify in detail the requirements on the mass of the disk and on the accretion rate posed on the 
models by realistic situations, we compare the SEDs produced by these models with the observed SEDs of a small 
sample of well-studied protostellar objects. We find that relatively modest disks - even lighter than the central 
star - can lead to an interesting fit to the infrared SED of the FU Orionis objects considered, while in the case 
of T Tauri stars the required parameters fall outside the range suggested as acceptable by the general theoretical 
and observational scenario. On the basis of the present results, we may conclude that the contribution of a self- 
gravitating disk is plausible in several cases (in particular, for FU Orionis objects) and that, in the standard 
irradiation dominated disk scenario, it would help softening the requirements encountered by Keplerian accretion 
models. 
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1. Introduction 

It now seems well established that pre-main-sequence 
stars, such as T Tauri and FU Orionis objects, are gen- 
erally surrounded by a disk of gas and dust. The tradi- 
tional evidence for such circumstellar disks comes from 
the excess luminosity found at long wavelengths, far- 
infrared (Rydgren, Strom & Strom 1976) and millimet- 
ric (Beckwith et al. 



199C), indicating that this radiation 



comes from a region of relatively low temperature, far from 
the central star. 

The first attempt at trying to explain such excess in- 
frared luminosity dates back to the 70's (Lynden-Bell & 
Pringle 1974), when it was argued that it might be due 
to the release of gravitational binding energy from an ac- 
cretion disk. In recent years, the general picture that has 
emerged is that, in most T Tauri systems, the mass accre- 
tion rate is too low (M ~ 1CU 8 — 1O _7 M0 /yr) to account 
for the observed luminosity and that the observed long 
wavelength excess is likely to result from the reprocess- 
ing of the star-light from the outer disk (Hartmann 1998Q . 
Therefore, the accretion disk of T Tauri stars is generally 



thought to be "passive". In contrast, during FU Orionis 
outbursts, the accretion rate appears to be relatively high, 
so that here the disk can dissipate sufficient amounts of 
energy to produce the observed luminosity. In the con- 
struction of dynamical models, a related issue is that of 
the total disk mass. Clues have indeed been gathered that 
the disk mass may be large in some protostellar systems 
(see Lay et al. 1994), especially in the earlier phases of the 
process of star formation (Yorke & Bodenheimer 1999, for 
example, conclude from hydrodynamical simulations that 
the disk mass at the beginning of the accretion phase is 
comparable to that of the central object) and in the con- 
text of high- mass star formation, where observations point 
to the presence of massive disks (Cesaroni et al. 



1994) 



Standard accretion disk theory (Shakura & Sunyaev 



1973) is based on the assumptions that the disk is optically 
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thick and geometrically thin and that the power viscously 
dissipated in the disk D = v a a(rQ,') 2 /2 is all radiated 
away. In the optically thick case, it can be shown that the 
spectral index n = dLog(^i 7, y )/dLogz/ in the infrared part 
of the spectrum is related to the exponent q in the surface 
temperature profile T s = Tpfro/r) 9 , so that n ~ 4 — 2/q 
(Adams, Lada & Shu 1988 ). For a Keplerian rotation curve 
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one has T s (r) oc r 3 / 4 , leading to L v = vF v oc ^ 4 / 3 , while 
for T Tauri stars one typically finds L v oc ^ 2 / 3 (Kenyon 



& Hartmann 1995). The difficulty becomes even more ev- 
ident for the so-called "flat-spectrum T Tauri" (such as 
T Tau itself and HL Tau), which show a spectral index 
n ~ 0. In turn, if the rotation curve is flat, the surface 
temperature profile decreases more slowly T s (r) oc r -1 / 2 . 
In this latter case, the spectral index would be n = 0, 
rese mbling that of flat-spectrum objects. Adams, Lada & 
Shu ( 1988 ) indeed proposed that the flat spectrum could 
arise from a disk with such a non-standard surface tem- 
perature profile, but argued that the observed disk mass 
in these systems (although difficult to measure) is proba- 
bly too low to produce substantial changes to the rotation 
curve and then suggested that the non-standard temper- 
ature profile should be ascribed to collective processes of 
angular momentum and energy transport associated with 
self-gravity induced instabilities. 

In the absence of a convincing physical justification 
for non-standard descriptions of accretion disks, the at- 
tention has been turned to other ways of explaining the 
far-infrared spectra of T Tauri stars. As noted above, it is 
now commonly thought that in most of these systems the 
accretion rate is low and that the heating of the disk is pri- 
marily due to irradiation from the central star. The effect 
of irradiation is enhanced at large radii by the disk flaring, 
when the thickness of the disk increases rapidly with ra- 
dius. Models of the vertical structure of these flared disks 



have been constructed (Kenyon & Hartmann 1987; Chiang 
& Goldreich 1997; but concerns about stability have been 



expressed by Dullcmond 2000). On the other hand, in the 



case of flat spectrum sources and of most FU Orionis ob- 
jects, flaring alone is ineffective and the far-infrared emis- 
sion has been attribut ed to an infalling envelo pe of dust 
(Kenyon & Hartmann |l99l| ; Calvet et al. |l994p . 

The hypothesis that the flat long wavelength spec- 
trum of protostellar objects arises from a self-gravitating 
disk with flat rotation curve actually remains attractive 
from a physical point of view, especially for its simplic- 
ity. Unfortunately, only scattered efforts have been made 
in the direction of constructing self-gravitating models of 
accretion disks. The general ideas mentioned above have 
been revisited recently along with a new viscosity pre- 
scription for self-gravitating disks (Duschl, Strittmatter 
& Biermann 200C), but little or no physical justification 
has been provided for the desired radial surface density 
profile. In view of this discussion, it seems natural to 
consider and to put to a test the accretion disk models 



that we have constructed recently (Bertin 1997; Bertin 



& Lodato 1999, hereafter BL99), which include the grav- 
itational fields of both the central star and the disk. 
Interestingly, these models show that the rotation curve 
of a self-regulated accretion disk is Keplerian in the inner 
disk and flat in the outer disk, corresponding to a surface 
density a oc r _1 (which, incidentally, has been observed 
m some cases: e.g., TW Hya, Wilner et al. |2000| , and HH 
30, Burrows et al. 1996| ), the transition taking place at a 
radius r s w (GM+/8)(GM /2a)~ 2 / 3 (using standard nota- 



tion; see Section g). Deviations from Keplerian rotation 
occur even for disk masses Mdisk ^5 which could be 
a plausible range even for some T Tauri disks. Based on 
these models of self-regulated accretion, it is now possible 
to quantify in detail what would be the parameter require- 
ments needed to fit the infrared spectra of protostellar 
objects. In this respect, our results will turn out to be 



quite different (see discussion in Sect. 3.5) from earlier es- 
timates ( e.g., s ee Kenyon & Hartmann 1987; Shu, Adams 



& Lizano 1987 ) that did not include the important ingredi- 
ent of self-regulation, which is likely to change significantly 
the energy balance equations in the se lf-gr avitating part 
of the disk (see also discussion in Sect. 2.2). 

In this paper we describe the spectral energy distribu- 
tion produced by a class of self-regulated, self-gravitating 
accretion disks and discuss its dependence on the various 
parameters involved and the differences from the non-self- 
gravitating case. In addition, we determine the physical 
parameters needed to fit the SEDs of realistic protostel- 
lar objects, by referring to a small sample of well-studied 
FU Orionis systems and T Tauri stars. A general theory 
of FU Orionis and T Tauri systems, or a comprehensive 
picture of the individual objects picked here to test the 
pure self-gravitating model against observed SEDs, is be- 
yond the goals of this paper. The models considered are 
extremely simple and do not address the issue of the phys- 
ical ingredients that determine the long-term evolution of 
the system (in particular, the physical ingredients respon- 
sible for the outburst phase of the FU Orionis systems): 
they are taken to match current "snapshots" of such evolv- 
ing systems, assuming that the steady-state equations are 
temporarily adequate to describe the observed situation. 

The quality of the fit to the available data for the long 
wavelength spectral energy distribution is very good. In 
the case of FU Orionis systems the required accretion rate 
is of the order of M w 10~ 5 — 10 _4 M Q /yr and the im- 
plied disk mass is of the order of one solar mass. The im- 
plied mass accretion rate for T Tauri stars is found to be 
M « 10" 7 -10" 6 A/ Q /yr. The choice of the outer radius of 
the disk determines the mass of the disk relative to that of 
the central star. The estimates, although uncertain, point 
to rather massive disks {Mdisk ~ 0.3M Q ). In the case of 
FU Orionis systems the mass accretion rates are in rea- 
sonable agreement with the currently accepted estimates, 
while the disk masses are higher than the commonly as- 
sumed disk masses in these systems. For T Tauri stars, 
instead, the required parameters tend to fall outside the 
range suggested as acceptable by the general theoretical 
and observational scenario. In any case, the contribution 
of self-gravity may soften the demands on models that at- 
tempt a description without the use of the disk self-gravity. 

The paper is organized as follows. In Section ^ we de- 
scribe the model of self-gravitating, self-regulated accre- 
tion disk that we adopt here; in Section || we describe the 
spectral energy distributions derived on the basis of our 
model and discuss their general properties; in Section |]we 
consider a sample of FU Orionis objects and fit their SED 
with the self-gravitating disk model, deriving the physical 
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properties of the disks inferred on the basis of our model, 
and discuss the results obtained; in Section |^ we point 
out the difficulties met when trying to adopt a similar ap- 
proach for the case of T Tauri stars; in Section ^| we draw 
our conclusions. 



2. Self-regulated accretion disks 

2.1. Disk model 

We refer to a self-gravitating accretion disk model that is 
only partially self- regulated (BL99), so that the relevant 
Q-profile is of the form: 



Q 



C S K 

TrGa 



1 



-9/8 



exp(-r/r Q ) 



(1) 



where c s is the effective thermal speed, k is the epicyclic 
frequency, and a is the surface density of the disk. In this 
way the stability parameter Q is much higher than unity 
in the inner disk (r < tq) and becomes constant (Q ~ 
Q w 1) in the outer disk. While in the context of galactic 
dynamics the scale tq is usually associated with the bulge 
size and, in the context of AGN, some efforts have been 
made in the description of the transition from a non-self- 
gravitating to a self-gravitating accretion disk (Bardou et 
al. 1998 ), no clear-cut prescription for the onset of self- 



regulation in protostellar disks is available. We have thus 
decided to take xq — rQ/r s (for the definition of r s , see 
Eq. (J2J) below) as a free parameter of our model. The scale 
r s is the scale that defines the transition from Keplerian 
to flat rotation curve: 



2GM* 



GM 
2a 



-2/3 



(2) 



Here M* is the mass of the central object, M the accre- 
tion rate, and a the viscosity parameter of the Shakura & 
Sunyaev ( 1973| ) prescription 



v a G — ac s h<j = —p;h(Q, Q/k); 
ttG 



(3) 



in the last equation we calculate the disk thickness as h = 
(c 2 s /TrGa)h, with: 



(4) 



4Q 2 [2fi 2 /K 2 - 1] 




The last expression gives a useful interpolation formula for 
the thickness of the disk in the case when the contributions 
of the star and of the disk are both taken into account in 
the vertical gravitational field (see BL99, Appendix). 

In the following, for simplicity we assume that the con- 
servation law for the angular momentum is written as: 



GMr 2 n + 2ac\hr 



dSl 







(5) 



and thus avoid the relatively complex behavior that is 
associated with a net transport of angular momentum 



(J 7^ 0; calculated by BL99). Here we may also recall 
that Popham et al. ( 1996| ) find better agreement with the 
observational data when their boundary layer models for 
FU Orionis systems are taken with a negligible value of J. 

Finally, the rotation curve is computed by solving the 
full Poisson equation, including the contribution of both 
the central star and the disk: 



, GM* 
V 2 = 



dr 



(G) 



where the gravitational field of the disk d^/dr is given 
in Eq. (4) of BL99. In practice, a useful way to express the 
rotation curve is by scaling V 2 to its value in the outer 
disk: 



2 \Q 



GM 
~2a 



2/3 



(7) 



with cj) 2 given in Eq. (13) of BL99. 

Equation (|J) gives an expression for the profile c s (r) in 
terms of HjV) and Q(r); the latter functions are available 
from Eq. (@) and Eq. ([!]). Similarly, Eq. (||) provides the 
profile of the quantity v a a . Therefore, we have an explicit 
expression, dependent on few parameters (such as a, M , 
Q), for the viscous dissipation rate D(r) = v a a(rVL') 2 /2. 
In the standard theory this would be sufficient to specify 
the surface temperature profile T s (r) associated with the 
optically thick emission of the disk. Here , as described in 
a separate article (Bertin & Lodato 2001, hereafter BL01; 
see also additional remarks in Sect. [2.2j below), it is ap- 
propriate to consider an additional heating term in the 
energy balance equation, so that the surface temperature 
profile is determined by: 



<jbTs = -v a a{rSl'Y 



9{Q) 



Q 2ttG 



(8) 



with g(Q) = (Q/Q) m , m a large number (we will take 
m = 20), and where ub is the Stefan-Boltzmann con- 
stant. The additional Jeans-related term ensures that the 
process of self-regulation properly occurs. Note that in 
Eq. (||) we neglect the contribution to the heating of the 
disk associated with the reprocessing of the light emit- 
ted by the central star, which is commonly thought to be 
the key heating source for T Tauri disks, especially if the 
outer disk is "flared"; we will discuss the effects of disk 



irradiation below (Sect. 3.4). In any case, our model of 
self-gravitating disk has negligible flaring (see BL99). 

As a final remark, we note that the assumption that 
the emission of the disk is optically thick is not conclusive, 
primarily because of the uncertainties in dust opacity. In 
this respect, our recent analysis of the energy budget in 
self- regulated disks (BL01), although it was focused on 
the different context of AGN, has pointed out that self- 
gravitating disks may be optically thin. Here, for simplic- 
ity and for a direct comparison with other studies, we 
consider models characterized by optically thick emission. 

The notation in the present paper is slightly different 
with respect to earlier articles (in relation to the use of 



4 



G. Lodato and G. Bertin: The SED of protostellar disks 



c s and v a ), in order to avoid confusion with the universal 
notation for some physical quantities (speed of light and 
frequency) . 

2.2. Self-gravity, viscosity, and the energy budget 

Where does the energy of the Jeans-related term in Eq. 
(||) come from? Why is the heating associated with insta- 
bilities not necessarily related to M (through v a u)l The 
reader is referred to a separate article (BL01) devoted to 
introducing and justifying our approach. Here we take this 
opportunity to add some additional comments that might 
help answer questions such as those posed above. 

Earlier attempts at incorporating the role of collec- 
tive effects associated with the disk self-gravity generally 
focused on the problem of angular momentum transport 
and thus on a modification of the viscosity prescription; 
these approaches (in part icular, see Lin & Pringle 1987, 



199C ; Bardou et al. 1998 ) often realized the need for a 



mechanism of self-regulation, locally sensitive to the value 
of Q. In turn, n umer ical simulations (e.g., s ee La ughlin 
& Bodenheimer 1994 , Laughlin & Rozyczka 1996 ) have 
demonstrated that spiral density waves can indeed trans- 
port angular momentum so as to act as an "effective vis- 
cosity" . In this respect, we agree that gravitational insta- 
bilities may have the welcome role, in the outer parts of ac- 
cretion disks, of contributing significantly to the effective 
viscosity of the disk, and hence directly to the accretion 
rate M. 

In the approach that we have followed (starting with 
Bertin 1997| ) we have argued that the process of self- 
regulation, which is expected to maintain the outer disk 
at values of Q of order unity, should result from a modifi- 
cation of the energy equations, so as to take into account 
the heating induced by gravitational instabilities. The pro- 
cesses involved are complex (see Sect. 3 in BL99 and Sect. 
3 in BL01). We should not expect that all the heating rate 
should be proportional to the angular momentum trans- 
port (as the ideal fluid equations with a viscous term gen- 
erally imply), especially since axisymmetric instabilities 
are known also to take place and gravitational heating is 
observed even in collisionlcss disks (e.g., see the simula- 
tions by Hohl 1971 and Hohl 1973). Therefore, we have 



thought it appropriate to separate the two issues, of the 
angular momentum and of the energy transport. 

In some extremely simple models (e.g., those consid- 
ered by N-body simulations of collisionless disks that bring 
out the effects of evolution induced by Jeans-related in- 
stabilities) one may argue, for example on the basis of 
the virial constraint, that if, globally speaking, energy 
is lost by the system, it must ultimately come from the 
gravitational potential well. Locally, the statement need 
not be true, given the facts that some energy may be re- 
distributed across the syste m thro ugh global instabilities 
(see remarks by Adams et al. 1989 ) and that gravitational 
energy has little local meaning. Thus, if we refer to Eq. 
(H), this local relation does neither contradict nor confirm 



the requirements of the global energy budget. (Note that, 
much in line with other descriptions of accretion disks, we 
do not make explicit use of a horizontal energy transport 
equation.) 

This general message, that the entire system cooper- 
ates globally to produce effects that show up differently 
in different parts of the disk, is often apparent from the 
results of a variety of numerical simulations (e.g., see 

200C| ). In par- 



Laughlin & Rozyczka 1996; Nelson et al 



ticular, several numerical studies have addressed directly 
the issue of the energy budget when dynamical instabil- 
itie s are involv ed (for some rece nt studies, see Pickett et 
al. |1998| , |2000t and Nelson et al. |l998| , |200CD . 

In addition, the disks that mediate the accretion pro- 
cess are actually not isolated: in fact, for the young stellar 
objects considered in this paper there is a central star with 
which the disk connects through a possibly hot and radia- 
tively inefficient inner boundary layer (for FU Orionis ob- 
jects, see Popham et al. 1996 ) and the disk is embedded 
in some environment from which the protostellar cloud 
has originated. These components, external to our set of 
model equations, represent a significant energy reservoir. 

In conclusion, we are not yet ready to produce an ex- 
plicit set of closed equations to describe the overall energy 
budget. In fact, to do so we should not only describe the 
global processes involved in the energy transport and the 
physical mechanisms active at the boundaries, but also 
face in detail the way we visualize an extremely complex 
system, in which the relevant heating and cooling terms 
may not be described by a simple one-fluid model (for 
example, for what concerns possible inelastic collisions 
between clumps). In this respect, it is interesting to see 
that in many numerical experiments relatively complex 
systems are generat ed, w ith transient interacting struc- 
tures (Pickett et al. 2000 ). On the other hand, the sim- 
pler heuristic approach that we have undertaken (in line 
with the spirit of the Shakura & Sunyaev prescription) 
can overcome this difficulty and provide insight into the 
properties of self-gravitating disks. 

While more work is certainly needed to set up a satis- 
factory global model, for the time being we proceed with 
this set of model equations that appear to possess most 
of the necessary ingredients required to describe a truly 
self-gravitating disk. 

3. General properties of the spectral energy 

distribution of active, self-gravitating accretion 
disks 

3.1. Spectral energy distribution for disk dominated 
cases 

From Eqs. (||), (||), and (||), it is useful to define a tem- 
perature scale To (with frequency isq — kT^jlixK) given 
by: 



n = 



M / 2 \ / M 



27rcr B a 2 \QJ 1 Mj ' 



(9) 
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with a dimensionless temperature T — T s /Tq and a di- 
mensionless frequency — 2ithv / kT$ = v /vq. 

For a disk inclined at an angle 9 with respect to the 
observer's line of sight, at a distance D , the luminosity 
is provided by integration of the Planck spectrum at the 
local surface temperature T s (r) between the inner radius 
ri n and the outer boundary radius r out of the disk: 

4,r/;i' ! 2;t/'(1/- 



vF„ = F 



cos 9 



n 2 



e 2irhv/kT s (r) _ I 



Using the dimensionless quantities defined above, the 
resulting spectrum can be written as: 

"* xdx 



L D 0" 



where: 



327r d rjfc>g 



COS( 



„6/T(x) 



240 cos6> 
^Q 2 



1 



M 



GM 
~2o7 



2/3 



(11) 



(12) 



Note that the scale Ljj of the disk luminosity does not 
depend on M*. 

The spectrum of the disk is then completely defined 
by two scale parameters Ljj and To and by four dimen- 
sionless quantities: the dimensionless inner and outer radii 
of the disk, x in = r ln /r s and x out = r out /r Sl the viscos- 
ity parameter a, and the dimensionless radius xq, defin- 
ing the inner boundary of the self-regulated region (see 
Eq. (0)). The two geometry independent scale parame- 
ters, Ljjj cos 9 and To can be traced back to the values of 
M and M*. 



Geometry 




cos 6 


= 1 


Accretion parameters 


M = lO" 5 M /yr 


M„ = 


1M Q 


Size of the disk 


Tin = 5Rq 


Tout — 


10AU 


Viscosity and self-regulation 


a = 10" 3 


XQ = 


= 0.4 



Table 1. Parameters of the reference model 



-1AU 
-5AU 

= 10AU (reference model) 
= S0AU 

= 10AU (Keplerian) 



° 32 




Figure 1. Shape of the spectral energy distribution for dif- 
ferent values of the outer radius r out . The quantity uF u is in 
arbitrary units. 



3.2. Effects generated by variation of the key 
parameters 

In this subsection we briefly describe the changes in the 
SED induced by changing some key parameters. We start 
from a specific reference model, for which the relevant pa- 
rameters are listed in Table |l|. For such reference model 
the scalelength r s turns out to be r s w 14AU, so that 
the scalelength beyond which the disk is self-regulated is 
r Q = x Q r s ~ 5.6AU. 

Figure [l] illustrates the effects of varying the outer ra- 
dius of the disk. As expected, increasing the outer ra- 
dius leads to a higher infrared excess. For comparison, 
Fig. [| also shows a model indicated by "Keplerian" . This 
model is computed with the same physical parameters as 
for the reference self-gravitating model (i.e., the same M, 
the same inner and outer radii, etc.), but as if self-gravity 
were turned off, i.e. by taking a strictly Keplerian rotation 
curve and by neglecting the additional term proportional 
to g(Q) in Eq. (||). The self-gravitating model clearly 
shows a much flatter SED in the far infrared. In con- 
trast, for r out — 1AU the two models (self-gravitating and 
Keplerian) would produce practically coincident SEDs. 
This is because for r out = 1AU the outer radius of the 
disk is smaller than tq, so that the disk has not become 
self- gravitating in this case. 



r m =2.5R 




12 13 14 15 16 



log(i//Hz) 

Figure 2. Shape of the spectral energy distribution for dif- 
ferent values of the inner radius fin- The quantity vF v is in 
arbitrary units. 

Figure ^| shows the effects of changing the inner ra- 
dius of the disk. Not surprisingly, this parameter influ- 
ences only the optical part of the spectrum. In Fig. |] we 
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x q =0.4 (reference model) 

x a -0.6 
x.=0.9 




log(i//Hz) 

Figure 3. Shape of the spectral energy distribution for dif- 
ferent values of the transition dimensionless radius xq. The 
quantity vF v is in arbitrary units. 



increase the value of xq. For self-gravity to give a sub- 
stantial contribution to the spectral energy distribution 
it is necessary that a significant portion of the disk be 
self-gravitating. For example, for xq = 0.6, which cor- 
responds to tq 8.5AU (recall that the outer radius is 
Tout — 10AU), the SED is still very close to that of the 
"Keplerian" case. 

3.3. Scaling relations with respect to a 

In the previous subsection we have described the changes 
in the resulting spectral energy distribution in response to 
changes in the physical parameters of the disk. Here, we 
will briefly discuss the constraints posed on the relevant 
parameters, when we consider the problem of fitting the 
SED of a given object by the self-gravitating accretion 
disk model. 

The highest temperature of the disk, which defines the 
typical frequency of the optical-near infrared "bump" of 
the SED, is determined, as in standard accretion disks, by 
the mass accretion rate, through the relation: 



3GMM+ 
8irr 3 ' 



(13) 



independent of the value of a. In a similar way, the typi- 
cal luminosity of the optical-near infrared part of the SED 
should be given by the standard non-self-gravitating ex- 
pression: 



GMM* 

27*7.77, 



(14) 



Therefore, if the SED of a given object (say, an FU Orionis 
object) is to be reproduced by a simple disk model, the 



observed typical wavelength and luminosity of the optical 
part of the SED (together with an estimate of M*), deter- 
mines the value of M and Ti n (the latter quantity turns 
out to be of the order of a few stellar radii; Hartmann & 



Kenyon 1996). On the other hand, in the self-gravitating 
case, the outer disk temperature, which defines the typical 
frequency of the long wavelength excess, is expected to be 
associated with a temperature scale T^ ut wTg/a (see Eq. 
(||)). In this way, the simultaneous constraints of fitting 
the optical and the far infrared parts of the SED define 
both a and M . 

Conversely, for objects such as T Tauri stars, for which 
the optical SED is dominated by the central star (see Sect. 
|3.4| below), the constraints set by Eqs. ( |l5| ) and ( |l4|) do not 
apply. In turn, under the hypothesis that the far infrared 
excess is due to the self-gravitating disk (which, however, 
may not be the case, because irradiation may often domi- 
nate), the temperature T^ ut rs T^/a is still constrained by 
the observations along with the luminosity scale of the far 
infrared spectrum, which is set by Ljj. In this case, when 
no further information on the accretion rate is given, for a 
given observed object, from Eq. (^), Eq. (^), and Eq. (12) 
we find the following scaling relations with respect to the 
viscosity parameter a: 



Mot a, (15) 

r s oc (M /a)~ 2/3 = const. (16) 

cos6> oc M^iM/a)- 2 ' 3 oc a -1 . (17) 

Indeed, in Section || we will see that, for two of the 
T Tauri objects considered, the exact value of a is not 
constrained solely by the observed far infrared excess, so 
that there is room for using different values of a. In the 
other cases shown, the disk gives a contribution also to 
the higher frequency part of the SED, and so there is little 
leverage on a, as in the disk dominated case. 

3.4. Inclusion of the central star and of disk irradiation 

During FU Orionis outbursts the luminosity of the disk 
is much higher than that of the central star, so this lat- 
ter contribution is generally negligible. In contrast, in the 
case of T Tauri objects the situation is reversed, so that 
the stellar luminosity is much higher than that of the ac- 
cretion disk. There can be however mixed situations in 
which both contributions should be taken into account. 
In this case, we should add to the right hand side of Eq. 
( |l0[ ) the contribution from the central star, which, for the 
limited purpose of this paper, can be taken to be that of 
a blackbody with given temperature: 

nRlinHu 4 1 

(18J 



D 2 



For convenience, we define a frequency — kT i ,/2nh. The 
resulting spectrum in terms of the dimensionless quanti- 



ties defined in Sect. 3.1 is then given by: 

4 -4 



9 15L* 
4^D 2 ^ = — j± 



-L D P 4 



xdx 
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where: 

L* = 4kH*<t b T*. (20) 

In addition to its direct contribution, the presence 
of the star affects the spectral energy distribution also 
through the reprocessing of the starlight by the disk. For 
a flat disk, the integrated effect of this contribution can 
be as high as L+/4. The effects of disk irradiation can be 
calculated by following Adams, Lada & Shu ( 1988 ), as an 
additional contribution to the heat balance equation, so 
that a term Fo (r) should be added to the right-hand side 
of Eq. (§), with: 



F D (r) = 



arcsm 




2a B T? (R* 



3tt 



(21) 



for a flat disk, and where the approximation is valid for 
r 3> R*. In terms of the dimensionless quantities defined 
previously, we find: 

4 / r, \ 3 

I n 

(22) 



F, 



D 



2 

3^ 



a 



<J B T* 

The first term on the right-hand-side of Eq. (@), associ- 
ated with viscous dissipation, gives a dimensionless con- 
tribution to the surface temperature profile ~ 1/x 3 , so 



that the parameter a defined in Eq. (22) measures the rel- 
ative importance of irradiation versus viscous dissipation. 
On the other hand, at large radii, when self-regulation be- 
comes important, the leading term that determines the 
surface temperature is the second, Jeans related term, on 
the right-hand-side of Eq. (||) . It can be easily seen that 

^ < 23 » 

where i?*, M*, and are the stellar radius, mass, and 
luminosity, respectively. As expected, the relative impor- 
tance of the stellar irradiation grows with increasing stel- 
lar luminosity and decreases with increasing mass accre- 
tion rate. 

In the following, we will generally neglect this irradi- 
ation term. We will check a posteriori the consistency of 
this assumption in the fitting procedure in the case of T 
Tauri stars, by calculating a on the basis of the fitted pa- 
rameters. 

In the case of disk flaring, it has been shown (Ruden 
& Pollack 1991) that the effect of irradiation is enhanced 
as: 



F D (r) 



3tt 



1 



h 



dln/i 



- 1 



(24) 



dlnr 

in the limit h, i?^ <C r. The term (dln/i/dlnr — 1) is deter- 
mined from the detailed vertical structure of the disk; its 
value is equal to 1/9 in the model by Kenyon & Hartmann 
( 1987 ) or to 2/7 in the model by Chiang & Goldreich 
( 1997 ). In our model (see BL99) at large radii the open- 
ing angle h/r does not depend on r and hence there is no 
flaring. 



3.5. Simple estimates of M and M^isk 

Here we will briefly discuss the simple arguments often 
provided (see Kenyon & Hartmann 1987 ) to estimate the 
values of the accretion rate M and of the disk mass 
Mdisk directly from the observations of T Tauri stars. 

In the case in which viscous dissipation is the only 
source of heating, the disk surface temperature can be ex- 
pressed approximately as: 



v B Tt{r) = -rft 



m dn Am 2 



4-7T dr 



4tt 



din ft 



dlnr 



GMM tot {r) 
47rr 3 



(25) 



where, in the last expression, M to t — + Mdisk- This 
expression, apart from numerical factors of order unity, 
is valid both in the Keplerian case (where M to t ~ -M* 
and ft 2 w GM+/r 3 ) and in the disk dominated "flat ro- 
tation curve" case, where Mdisk{r) oc r and ft oc 1/r. 
For a disk with a transition from an inner Keplerian re- 
gion to an outer self-gravitating region with flat rotation 
curve, the mass of the disk in its outer part is given by 
Mdisk(r)/M+ fts r/r s , where r s is the transition radius (see 
Eq. (|). 

If we try to follow the argument by Kenyon & 



Hartmann (1987), for the case of T Tauri stars, in or- 
der for the outer disk to produce the required infrared 
excess at wavelengths of the order of 25 — 100/im, its tem- 
perature should be of the order of 100K. Therefore, the 
required product MM to t{r ut) can be estimated by fixing 
T s (r out ) = 100K (with r out w 20AU; but note that a value 
of r out ~ 10AU would be more appropriate if we refer to 
the models described in Sect. || below). With this proce- 
dure one gets_MM tot /(10~ 7 M|yr- 1 ) » 10 3 . Kenyon & 



Hartmann ( 1987 ) then noted that both alternatives left 
by an interpretation in terms of an active disk would be 
highly implausible: (i) the accretion rate is high, at least 
in the outer disk, while the disk mass remains very small, 
thus leading to far too high values of the accretion rates 
(M > 10~ 5 M®/yr), or («) the disk mass at r out is very 
high, of the order of 1000M Q , while the accretion rate re- 
mains small (M < 10~ 7 M Q /yr). In this latter case, to 
have such a high disk mass, the transition radius to the 
self-gravitating part of the disk should be very small. In 
fact, according to Kenyon & Hartmann (1987) the transi- 
tion would occur at r — i?*. 

In reality, in a truly self-gravitating case, viscous dissi- 
pation is not the only source of heating (see discussion in 
Sect. 2.2 and BL01), and so Eq. ( p5| ) cannot be used. In 
addition, when a detailed model of self-gravitating accre- 
tion disk of the kind adopted in this paper is considered, 
on sees that the transition radius is not an independent 
parameter. In particular, from Eq. (Q), when M* = O.6M 
and M = 10 -7 M Q /yr, we see that r s w 40AU, so that the 
implied disk mass need not to be so high. On the other 
hand, when self-gravity effects are fully incorporated, the 
surface temperature is going to be higher than expected 
from the simple estimates associated with Eq. (P5|). In 
fact, from Eq. (||), especially when a is small, it is clear 
that even with mass accretion rates lower than those in- 
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ferred from arguments based on Eq. (|2j^) for the "small- 
disk-mass" active disk, an active disk could produce the 
required luminosity. 

In Section || this issue will be discussed further, in view 
of the parameters determined by fitting concrete examples 
of spectral energy distributions of T Tauri stars. Here we 
would like to anticipate that, even though the results, as 
already suggested in this subsection, show that the re- 
quired mass accretion rates need not be as high as esti- 
mated earlier, for the majority of T Tauri stars, the picture 
of a pure, self-gravitating, active disk remains unsatisfac- 
tory. 



4. A possible fit to the spectral energy 
distribution of some FU Orionis systems 

In this Section we will compare the model described in 
the previous Sections with the observations by fitting the 
spectral energy distribution of some FU Orionis objects. 
In these systems there is already evidence that the accre- 
tion disk is active, i.e. that most of the emission from these 
objects is probably due to the disk itself, which outshines 
the central star. There are already clues that disk self- 
gravity should play a major role in these cases; in particu- 
lar, most outburs t mod els predict the existen ce of massive 
disks (Bell et al. |l995[ Hartmann & Kenyon |l996[) . 



4.1. A sample of FU Orionis systems 



star 


M*/M 


A v 


D 


FU Ori 


1 


2 


550 pc 


V1515 Cyg 


1 


2.8 


1000 pc 



Table 2. Stellar parameters assumed for the FU Orionis sam- 
ple stars. 



Currently, there is a handful of known FU Orionis ob- 
jects, among which the three best studied systems are FU 
Ori, V1515 Cyg, and V1057 Cyg. The last object shows a 
relatively rapid decline in B magnitude after the outburst, 
with a decay timescale of a few years, about ten times 
faster than FU Ori, and thus special care must be taken 
when comparing observations made at different epochs. 
For this reason we have not included V1057 Cyg in our 
sample of FU Orionis systems. In addition to FU Ori and 
V1515 Cyg, we have also studied V1735 Cyg, which is an 
"extreme" case, with a huge infrared excess, so that the 
infrared luminosity is apparently larger than the optical 
one. Actually, this object is heavily obscured, so that it is 
very difficult to derive the extinction correction with con- 
fidence and hence to extract for it the intrinsic spectral 
energy distribution. Therefore, although we initially kept 
this object in our sample, so as to have an idea of how 
large the disk mass should be when the infrared excess is 



so large, we do not record here the results of the fit for 
this case. 

The parameters to be determined by the fit are: the 
characteristic frequency v>q, the dimensionless inner and 
outer radii of the disk Xi n and x ou t, the inclination angle 9, 
and the viscosity parameter a. We have derived the best- 
fit parameters for different choices of xq taken to vary in 
the interval [0.4, 0.9]. To derive the physical parameters of 
the disk (such as the accretion rate and the disk mass), 
we have to specify the central star mass (see Table |). In 
the case of VI 735 Cyg, we have made a fit (not shown 
in detail here) by assuming Ay = 0, M* = 2M Q , on 
the high side, given the unusual infrared luminosity of 
the object (see also Sect. 3T). In any case, in contrast 



with the models that are generally considered, we do not 
include here irradiation by the inner on the outer disk (our 
outer disk has no flaring and hence this contribution is 
negligible) and the possible contribution of an additional 
inf ailing envelope. 

The optical and near-infrared luminosities are those 
collected by Kenyon h Hartman n (|l99l[ ). The IRAS fluxes 
are taken from Weaver & Jones (1992). First the data have 
been dereddened using the procedure of Cardelli, Clayton 



& Mathis (1989) and then the spectral energy distribu- 
tions described in the previous Section have been fitted to 
the data. 

The best-fit parameters are shown in Table || and the 
derived physical parameters of the disk are shown in Table 
||. Figure |J displays the spectral data of the stars together 
with our best-fit SEDs (solid lines) and the corresponding 
SEDs for non-self-gravitating, Keplerian disks (thin lines), 



with surface temperature T s 



-3/4 



In Fig. H we plot, 



as an example, the surface temperature profile and the 
rotation curve for the best-fit model to FU Ori (in the 
case xq = 0.4) together with the corresponding quantities 
for the Keplerian model, showing that the outer parts of 
the disk in the self-gravitating case are indeed hotter and 
rotate faster than the corresponding Keplerian disk. 

The mass accretion rates are of the order of 10~ 5 — 
10 _4 MQ/yr. The value of a determined by the fit is in 
the range 10~ 3 — 10~ 2 , consistent with some theoretical 
expectations. Typical disk radii are of the order of 50AU. 
Disk masses derived from the fit are of the order of one 
solar mass. Some of these conclusions would also apply 
to V1735 Cyg, except that in this case, as a result of its 
very high infrared excess, the estimated disk mass would 
be around 20M Q . 

4.2. Discussion 

The infrared spectral energy distribution of FU Orionis 
systems has been described in terms of active accretion 



disks b y Kenyon et al. (1988) and by Kenyon & Hartmann 
( 1991 ). These authors find good agreement between the 
accretion disk model and the observed spectra for wave- 
lengths smaller than 10/im. The excess emission at larger 
wavelengths is attributed by Kenyon & Hartmann ( 1991 ) 
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FU Ori 


V1515 Cyg 


Log(v / Hz) 


[12.00-12.09] 


[11.94-12.12] 


cos 6 


[0.65-0.72] 


[0.94-0.88] 


a 


[3.4 10" 2 -1.61 10 -2 ] 


[2.17 10" 3 - 1.09 10~ 3 ] 


X in 


[0.0013-0.0018] 


[9 10 _4 -0.0016] 


Xout 


[1.52-2.51] 


[1.04-1.96] 



Table 3. Best-fit parameters for two FU Orionis objects. The values in square brackets refer to different choices of the free 
parameter xq ranging in the interval [0.4 — 0.9]. 





FU Ori 


V1515 Cyg 


M/(10- 5 M Q /yr) 


[8.2-6.6] 


[1.1-1.2] 


M dlsk /M e 


[1.7-2] 


[0.91-1.52] 


r 3 /AU 


[38-26] 


[23-14] 


9 


[49° _ 44°] 


[20° - 28°] 



Table 4. Derived parameters for the FU Orionis objects. The brackets refer to the same choice of xq as in Table 



Wi 35 

u 

0) 



O 
O 



1 1 1 1 1 1 1 1 1 1 

FU On 








A 


V 




\ 


; \ h / 





Log(i//Hz) 




Log(iv/Hz) 



Figure 4. Best fit to the spectral energy distributions for the FU Orionis sample (solid curves) compared to the spectra 
r esult ing from Keplerian disks wit h the same parameters (M, cosf?, Xi n , x out , a; thin curves). Data from Kenyon & Hartmann 
(p.99l|) and from Weaver & Jones (Il99§). 



to a flattened envelope which reflects the luminosity of the 
disk. The hypothesis that this excess comes from a flared, 



outer disk is discarded by Kenyon & Hartmann (1991), 
because it would have required unreasonable levels of flar- 
ing. 

Our estimates of the mass accretion rates are very close 
to those by Kenyon & Hartmann (1991), because these 



numbers are determined mainly by the optical-near in- 
frared part of the spectral energy distribution (see also 
Sect. |jO|); indeed, as can be seen from Fig. [|, the spec- 
trum for A < 10/xm is not affected by self-gravity related 
effects. On the other hand, in contra st with the interpre- 
tation of Kenyon & Hartmann ( 1991 ), we attribute in this 
paper the far infrared emission to the self-gravitating part 
of the disk. The mass thus required to produce the ob- 
served excess is of the order of one solar mass, which is 



rather high, but in line with the indications of current out- 
burst models (Bell et al. 1995 ). Unfortunately, no separate 
measurements of the disk mass in FU Orionis systems is 
available to our knowledge. 

The interpretation in terms of the presence of an in- 
falling envelope responsible for the far infrared excess is 
actually appealing because it can also offer a clue to the 
peculiar time evolution of V1057 Cyg, which shows a cor- 
relation between the decay timescale in the optical and 
that at longer wavelengths (lO^tm). Here we do not ad- 
dress the time evolution of our model, and indeed we have 
preferred to leave V1057 Cyg out of our sample. Still, it 
would be interesting to consider these issues in the con- 
text of self-gravitat ing d isks (see also the outburst model 
of Armitage et al. 2001 , which includes some effects re- 
lated to the disk self-gravity). The far infrared SED in 
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E- 



FU On 



self-gravitating 
Keplerian 



FU On 



self-gravitating 
Keplerian 



r/AU 



r/AU 



Figure 5. Surface temperature profile (left) and rotation curve (right) for the best-fit model to FU Ori (with xq = 0.4), 
compared with the corresponding Keplerian model. 



our models is sensitive to the relative positions of the disk 
outer radius r ou t and the radius rq at which the disk be- 
comes self-regulated (see Sect. [T^). A detailed description 
of the role of self-gravity in the context of time-dependent 
outburst models, which might result from the interplay 
between tq and r out as a function of time, is, however, 
beyond the goals of this paper. 



5. Difficulties of the self-gravitating scenario in 
the context of T Tauri stars 

In the context of T Tauri stars the situation seems to be 
less favorable for disk self-gravity to play a major role, be- 
cause the infrared excess is likely to result in general from 
irradiation of the disk from the central star. Many points 
are often mentioned in favor of this interpretation, among 
which we may recall the explanation of the observed sil- 



estimated. Second, as pointed out in Sect. 3.5 and as the 



icate emission features at lO/im (Calvet et al. 1992) and 
other HST observations (e.g., see Stapelfeldt et al. 1998| ). 
On the other hand, initial efforts at explaining the SEDs 
of T Tauri stars with active accretion disks led to esti- 
mated mass accretion rates so high that the corresponding 
optical-infrared emission of the disk should outshine the 
central star, which is not observed. In addition, the pecu- 
liar "flat" shape of the far infrared SED, when explained in 
terms of a non-Keplerian disk, led to un reasonably high 
disk m asses (Kenyon & Hartmann 1 987j ; Shu, Adams & 
Lizano [l987| ). 

The model described in this paper presents significant, 
qualitative and quantitative, differences with respect to 
the non-Keplerian models considered in previous investi- 
gations. First, the additional heating term in Eq. (|J) gives 
an important contribution to the heating of the disk in its 
self-gravitating part, so that the corresponding mass ac- 
cretion rates needed to account for the infrared luminosity 
in an active disk scenario may be lower than previously 



fit to FU Orionis objects in the previous Section demon- 
strates, the disk masses needed to give a significant contri- 
bution to the SED in the self-gravitating scenario are not 
as high as previously argued. Therefore, we think it ap- 
propriate to reassess the issue of the required parameters 
within the self-gravitating active disk scenario. 

Furthermore, in some cases (see Gullbring et al. 2000; 
Testi et al. 2001), such as in some intermediate mass 



Young Stellar Objects, the disk mass may be relatively 
high, so that the self-gravity of the disk should play some 
role, even without abandoning the general fra mew ork of 
irradiated disks. In addition, as noted in Sect. 3J5, when 
M* 0.6M Q and M w lO~ 7 Af /yr (which is not unrea- 
sonable at least for some objects), the implied values of 
r s readily show that the contribution of disk self-gravity 
should be incorporated. 

For these reasons we now proceed to fit the SEDs of 
some T Tauri stars with the model presented in this pa- 
per. The primary goal is to check in detail the parameter 
requirements for the active, self-gravitating disk scenario, 
so as to confirm or to correct previous estimates, even 
though this model is going to be relevant only for a small 
set of less typical cases. 

5.1. A sample of T Tauri stars 

We have considered two typical T Tauri stars (BP Tau 
and DE Tau) in the Taurus- Auriga molecular cloud (at 
a distance Do = 140 pc) for which the stellar masses are 
available from Gullbring et al. (1998) and the stellar tem- 
peratures and visual extinction from Kenyon & Hartmann 
( 1995 ). In addition, to check the parameter requirements 
on our models set by less typical cases, we have consid- 
ered two "continuum" T Tauri stars, DR Tau and DG 
Tau, in which the high amount of veiling in the optical 
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suggests the existence of high accretion rates (Gullbring 
et al. 2000). DG Tau is also one of the "flat spectrum" 
T Tauri stars considered by Adams et al. ( [1988 ) . The 
stellar parameters for the last two cases have been taken 
from Gullbring et al. (2000), with stellar temperature 
Ti, = 4000K and assumed stellar mass of 0.5M©. 

To fit this set of T Tauri stars, we add to the disk emis- 
sion also the direct contribution of the star, as described in 
Sect. iA. Nevertheless, we still neglect the effect of irradia- 
tion. The consistency of this assumption will be checked a 
posteriori by estimating the parameter a from the derived 
fit parameters. This, of course, translates into an actual 
overestimate of the mass accretion rate, as derived by the 
fit. The luminosity of the star L+ is added to the set of 
parameters to be determined by the fit. 

The optical and near-infrared lumino sities for all stars 
are available from Kenyon & Hartmann (1995). The IRAS 
fluxes for all stars are taken from Weaver & Jones ( 1992| ). 
The dereddening procedure was the same as adopted for 
the FU Orionis objects in the previous Section. For BP 
Tau and DE Tau, we found that the value of the viscosity 
parameter a does not affect significantly the quality of 
the fit, so it has been excluded from the fit parameters, 
and the fit has been performed, as separate examples, for 
two different values: a = 10~ 3 and a = 10~ 4 (see also 



comment at the end of Sect. 3.3). 

The best-fit parameters (for different values of the pa- 
rameter xq, taken to vary in t he in terval [0.4,0.9]) are 
shown in the Appendix in Table A.l and the derived pa- 
rameters are reported in Table A2. The fitted SEDs are 
shown in Fig. |g. As can be seen, the fitted solid curves 
reproduce the available infrared data well. The surface 
temperature profile and the rotation curve for DR Tau 
(in the case where xq = 0.4) are shown in Fig. 0. The re- 
sulting mass accretion rates for BP Tau and DE Tau are 
of the order of 10~ 7 MQ/yr for a = 10~ 4 and larger by a 
factor w 10 for a — 0.001, as a ntic ipated from the scaling 
relations pointed out in Sect. |3.3[ In these two cases, it 
can be seen that increasing the value of xq (with a being 
fixed) leads to higher disk masses and accretion rates, and 
consequently to lower values of r s (which is proportional 
to M -2 / 3 ) and a (oc M _1 ). Typical disk radii r out are of 
the order of 10 — 50 AU. The masses required to fit the 
data are of the order of some fraction of M & . 

For DR Tau we obtain M « 5 10~ 7 M G /yr, while for 
DG Tau we have M » \Q- & M Q /yr. The best fit value of a 
is of the order of 10~ 4 and slowly decreases with increasing 
xq. The required values of the disk masses are a fraction 
of Mq also in these cases. 

5.2. Discussion 

For T Tauri stars, the infrared spectral energy distribution 
is generally not considered to be a good tracer of the mass 
accretion rate, given the fact that a significant fraction of 
the infrared luminosity is likely to be due to reprocessing 
of the starlight rather than to accretion power. Mass ac- 



cretion rates are better derived from the hot continuum 
expected to be produced from the accreting matter when 
it hits the stellar surface (Gullbring et al. 1998[ ). In any 
case, great uncertainties remain in the derived accretion 
rates (especially when the hot accretion continuum is very 
high) mainly due to uncertainties in the extinction correc- 
tions. The three most extensive samples of T Tauri stars 
with accretion rates determined in this way (Valenti et al. 



1993; Hartigan et al. 1995; Gullbring et al. 1998) report 



values that differ for the same star by as m uch as a factor 
of ten. On the averag e, the Gullbring et al. (|l998|) and the 
Valenti et al. ( 1993 ) estimates are consistent with each 
othe r, alt hough with large scatter, while the Hartigan et 
al. ( 1995 ) estimates are systematically biased to higher 
values. 



Looking at the values in Table A. 2 for the cases of BP 
Tau and DE Tau, it is clear that the pure self-gravitating 
disk scenario is unsatisfactory. For a = 10 -3 the values 
of the required accretion rates for these objects are of the 
order of lO _6 M0/yr, two orders of magnitude larger than 
the values typically quoted (Gullbring et al. 1998). This is 
not surprising, because it is known that an interpretation 
in terms of active disks requires high mass accretion rates. 
These high accretion rates also lead to unreasonably small 
inclination angles for the disks. One might imagine that 
photo-evaporation of the disk due to UV radiation field 
of the central star could reduce the mass accretion rate 



in the inner disk (Richling & Yorkc 1997), thus allowing 
for discrepancies between the M derived from the infrared 
SED (which is most sensitive to the conditions in the outer 
disk) and that obtained from the hot continuum. However, 
this process requires the presence of a strong UV source 
that is not available for the T Tauri stars that we consider 
here. On the other hand, as suggested at the beginning of 
this Section, the required disk mass is not as high as pre- 
viously argued, because it turns out to be just a fraction 
of the stellar mass. 

In the case where a smaller value for a is assumed, the 
required mass accretion rates, in line with the scaling re- 
lation of Sect. 3.3, are lo wer, b ut still closer to the higher 
values of Hartigan et al. ( 1995 ) than to the mutually con- 
sistent values of Valenti et al. (1993) and of Gullbring et 
al. (1998). In this case, however, the lower accretion lu- 
minosity is accompanied by an increasing importance of 
the irradiation term, quantified by the higher values of a, 
thus leading to a mixed case where both contributions (of 
irradiation from the central star and of internal heating of 
the disk) are important. Finally, note that the luminosity 
of the star L* is smaller than usually estimated, mainly 
because we attribute part of the near-infrared luminosity 
to the disk, rather than to the star. 

The interplay between standard models and the self- 
gravitating scenario might be more interesting for the 
"continuum" stars, for which the best fit accretion rates 
from the model studied in this paper turn out to be higher, 
but not dramatically, than what inferred fro m obs erva- 
tions at shorter wavelengths (Gullbring et al. 2000 ). We 
should reiterate that, by our modeling procedure, we ac- 
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Figure 6. Best fit to the spectral energy distributions of our sample stars (solid curve), compared to the spectra resulting from 
Keplerian disks with the same parameters (M, co sff, X j n , x out ; thin curves). For al l cases the displayed fit refers to xq = 0.4. 
The spectral data are from Kenyon & Hartmann (1995) and Weaver & Jones (1992). 



tually overestimate the mass accretion rate because the 
contribution of irradiation is neglected. 

Figure ^ shows a significant discrepancy between the 
model curves and the data at short wavelengths, particu- 
larly for DR Tau and DG Tau. Here we are generally not 
concerned about the short wavelength data, because the 
optical and ultraviolet emission and its variability are well 
interpreted in terms of magnetospheric accretion (e.g., see 
Kenyon et al. |1994[). 



6. Conclusions 

In this paper we have described a model for self- 
gravitating accretion disks in the context of protostellar 
disks. This model might in principle be relevant for a num- 
ber of cases, because the mass of the protostellar disk is 
likely to be high in the earlier phases of the accretion pro- 



cess. In addition, massive disks are sometimes observed in 
the context of intermediate mass star formation. 

We have described the properties of the spectral en- 
ergy distribution of these self-gravitating disks and their 
dependence on the various parameters involved. In addi- 
tion, we have checked the parameter requirements for self- 
gravitating disks to fit the SEDs of some observed cases. 
In particular, we have considered the cases of FU Orionis 
objects and of some T Tauri stars. 

For the FU Orionis objects we find that the infrared 
excess at wavelengths larger than A = 10/im can be de- 
scribed in terms of the presence of a self-gravitating disk. 
Our models are otherwise consistent (for the implied mass 
accretion rate, for ex ample ) with previous models (see 
Kenyon and Hartmann 1991), which, in contrast, attribute 
the far infrared excess to an inf ailing envelope. 
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We have also compared the spectral energy distribu- 
tion of our models with the SEDs of some T Tauri stars, 
even though in general these should not be considered for 
a realistic application of our models, in order to check the 
requirements on an interpretation based on an active disk 
when self- gravity is taken into account. The required disk 
masses appear to be much lower than estimated previ- 
ously. On the other hand, especially when the viscosity 
parameter is taken to be high, we do find unreasonably 
high required mass accretion rates. For typical T Tauri 
stars, our very simple models have an evident limitation, 
because we have not considered the important role of disk 
irradiation. 

In general, we would like to see this paper as an in- 
teresting starting point for the inclusion of self-gravity re- 
lated effects, even in the standard framework of the irra- 
diated disk model. 
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Appendix A: Fit parameters obtained from the 
SEDs of some T Tauri systems 

In this Appendix we report the detailed values of the fit 
parameters and of the physical quantities derived from 
the fit of some T Tauri stars presented in Section 0. Table 
|A.l| shows the fit parameters (for BP Tau and DE Tau 
the parameters are s hown for two adopted values of a; see 
Section ||) and Table AJ2 displays the derived parameters. 
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BP Tau DE Tau BP Tau DE Tau DR Tau DG Tau 



Log(z/ / Hz) 


[11.86-12.11] 


[11.55-11.81] 


[11.77-12.03] 


[11.43-11.68] 


[11.60-11.81] 


[11.68-11.87] 


Log(L*/ (erg/s)) 


33.298 


33.231 


33.298 


33.231 


33.562 


33.499 


Log(L D /(erg/s)) 


[29.80-30.14] 


[29.42-29.77] 


[29.67-30.02] 


[29.26-29.60] 


[30.34-30.63] 


[30.62-30.89] 




[0.004-0.009] 


[8 10~ 4 -0.002] 


[0.003-0.007] 


[7 10~ 4 -0.001] 


[0.003-0.005] 


[0.004-0.007] 


Xout 


[0.60-1.32] 


[0.75-1.63] 


[0.60-1.06] 

(io- 4 ) 


[0.56-1.25] 
(IO" 4 ) 


[0.54-1.09] 


[0.81-1.56] 


a 


(io- 3 ) 


(10~ 3 ) 


[2.2 10~ 4 -0.8 10" 4 ] 


[4.7 10~ 4 -1.8 10" 4 ] 



Table A.l. Best- fit parameters for the T Tauri sample. The values in square brackets refer to different choices of the free 
parameter xq ranging in the interval [0.4-0.9]. The values of a in parentheses are assumed (see discussion in Sect. 0). 



a = 10 -3 a = 10~ 4 





BP Tau 


DE Tau 


BP Tau 


DE Tau 


DR Tau 


DG Tau 


M/(10- 6 M Q /yr) 


[3-7] 


[0.81-1.7] 


[0.5-0.7] 


[0.12-0.26] 


[0.5-0.5] 


[1-1] 


M d i sk /MQ 


[0.28-0.44] 


[0.19-0.31] 


[0.28-0.3] 


[0.13-0.21] 


[0.25-0.33] 


[0.41-0.5] 


r 3 /AU 


[15-9] 


[20-12] 


[11-7] 


[16-10] 


[20-10] 


[19-11] 


L*/L@ 
a/W' 1 


0.51 


0.44 


0.51 


0.44 


0.95 


0.82 


[0.06-0.02] 


[0.3-0.2] 


[0.4-0.1] 


[3.5-1.5] 


[1-1] 


[0.4-0.4] 


cos 8 


[0.03-0.02] 


[0.14-0.08] 


[0.12-0.07] 


[0.53-0.31] 


1 


1 



Table A. 2. Derived parameters for the disks of the T Tauri sample stars. 
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